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Abstract. We consider the Riesz and Hardy-Littlewood wave i.e. a "critical 
function" whose behaviour is concerned with the possible truth of the Riemann 
Hypothesis (RH). The function is studied numerically for the case a = ^ and 
P = 4 in some range of the critical strip, using Maple 10. 

In the experiments, N = 2000 is the maximum argument used in the Mobius 
function appearing in cj, i.e. the coefficients of Baez-Duarte, in the repre- 
sentation of the inverse of the Zeta function by means of the Pochammer's 
polynomials. 

The numerical results give some evidence that the critical function is bounded 
for 9t(s) > 5 and such an "evidence" is stronger in the region SH(s) > j where 
the wave seems to decay slowly. This give further support in favour of the 
absence of zeros of the Riemann Zeta function in some regions of the critical 
strip (£K(s) > v) and a (weaker) support in the direction to believe that the 
RH may be true (<H(s) > §). 



The starting point of this note is the representation of the reciprocal of the 
Riemann Zeta function by means of the Pochammer's polynomials Pk(z) (where z 
is a complex variable), whose coefficients Cfc have been introduced by Baez-Duarte 
for the Riesz case (a = (3 = 2), For the study of the coefficients Cfc, some recent 
analytical as well as numerical results have been obtained [2, 3, 4, 5, 6, 7, 8, 9]. For 
a rigorous treatment of the Miintz formula to the finding of new zero free regions 
of the Riemann Zeta function, the reader may consult the work of Albeverio and 
Cebulla [1]. 

Using the Baez-Duarte approach, the representation of may be obtained for 
a family of a two parameter Pochammer's polynomials (parameters a and (3 [4]) 
and reads: 
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1. Introduction 
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(1-3) c k (a,P) : = E ^(1-^)* 



fe 

1 

The expression for Cfe we will use in our computation is given by: 

N 



(1.4) ^) = E$^> 

n— 1 

which for large A; is a correct formula for the Ck given above (see Appendix) . From a 
theorem of Baez-Duarte [2, 3], an important inequality concerning the Pochammer's 
polynomials is given by: 

(1.5) \P k (z)\ <Ak-^ z \ 

The inequality, when applied to our family of Pochammer's polynomials gives: 



(1.6) 



From this it follows [2, 3, 4] that the RH will be true, i.e. that in the 
representation above will be different from infinity (no zero of ((s) for 9i(s) > p) 

if Ckk 2 ^ < C. For the numerical study it is convenient to introduce the variable 
x = log(fc), in term of which we define the critical function corresponding to a and 
(3. This is given by: 

2000 , \ 

M») --4 



(1.7) iP(x;a,f3,p):=e^ x J2 



n=l 

2000 is the maximum argument N used in these experiments, which for the 
special case we treat (a = and /3 = 4) ip will be calculated up to x — 30 (this 
corresponds to k = e 30 = 1.06 8 65 x 10 13 ). 

Before we present the results of our numerical experiments for various values of 
p (for p — l)|)f)f)5)f)^) it is important to give the explicit expression of the 
contribution of the non trivial (ip n t) an d a l so of the trivial zeros (ipt) to the critical 
function defined above for the general case a and j3, following the expression given 
by Baez-Duarte for the case a = /3 = 2 [2]. For the non trivial zeros, in the variable 
x = log(fc), at p, it is given by: 

(1-8) il> nt (x;a,(3,p) = CW^ ' 

where z is any nontrivial zero of C( s )- I n our experiments we will limit to the con- 
tribution of the first two lower zeros given experimentally by z\ = \ + il4. 134725... 
and Z2 = 5 + i21. 022040... and the complex conjugate of them. The corresponding 
contribution will be denoted by r\{x) (from z\ and zi) and r2(x) (from Z2 and 22)- 
The contribution of the trivial zeros z — — 2n for every integers n, to the critical 
function is given by: 



j 00 c -^ £ XY^ a+2n ^ 



(U)) C'(-2 
where a summation until N = 20 will be sufficient. 
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So, in our calculations we will set a = ^ and (3 = 4 in the above formulas, for 
any value of p we shall consider. The contribution ip t for p will be indicated with 
g p {x). Below we present the results of our numerical experiments performed using 
Maple 10, where as anticipated the maximum argument in the Mobius function 
present in the definition of the critical function (essentially Ck), is N = 2000. The 
fluctuations errors around 2000 will be specified in the Appendix. 

2. Numerical experiments 

In Fig. 1 we give the plot of the two functions ip(x; ^,4, ^) — r\{x) — r 2 {x) and 
gi/2{x) up to x = 30 which shows a very good agreement. Notice that we have 
taken into account only the contribution of the first two nontrivial zeros in the 
Baez-Duarte asymptotic formula for the c k . For the Riesz case (a = f3 = 2), the 
contribution of the trivial zeros to the c k have been treated by Maslanka using the 
Rice's integrals [8]. 




Figure 1. Plot of ip(x) — r\{x) — r%(x) [red] and 171/2(2;) [green] 
up to x — 30 

In Fig. 2 we present the the plot of the two functions tp(x; 4, |) — gi/2(x) and 
r\ (x) + ri (x) up to x — 30 which shows not only a good agreement but also the 
oscillatory behaviour of the contribution of the first two nontrivial zeros. 

In the next Fig. 3 we present the plots of some critical functions (ip p ) corre- 
sponding to different values of p using (1.7) and this without any comparison with 
the Baez-Duarte asymptotic expansion considered above. It is to be noted that all 
functions ijjp has the same zeros and we observe that there is a well marked evidence 
that for p> \ increasing to 1 the amplitudes decay while for p < \ the amplitudes 
grow. These functions have been indicated with tpi, ip7/8, ^3/4 > V's/Sj ^1/2 > ^3/8' 
^3/10 respectively. 
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Figure 2. Plot of ijj{x) — gi/2(x) [red] and r 1 (x) + r 2 (x) [green] 
up to x — 30 

It should be said that ip^/s and 1P3/10, we have considered, have no relation with 
the representation of which is valid only for SH(s) > \. The two functions help 
only to visualize that ipi/ 2 ^ s tne borderline for the critical functions decaying for 
$H(s) > | as suggested by our numerical experiments up to x = 30. It should also 
be added that from the duality relation (Riemann symmetry of the Zeta function) , 
given by: 

/ \ i ...rfe) i 

it follows that the right hand side of (2.1) ensures a representation of via the 

Pochammer's polynomials in the region < 9t(s) < \- 

In Fig. 4 we present the result for a special case where we allow a slower decrease 
in the critical function (see addendum in the exponent of the critical function), 
which is the same as to say that we ask only for a slower decay of Cfc, at p = \ i.e. 
of the type Ck — Alo ^ k ' > for the case considered. This is not the same as to ask 

fc4 

that RH is true or that RH is true with nontrivial zeros which are simple [3]. It is 
a case in between the two. 

In this case the critical function (indicated with "01/2+) is explicitly given by: 



(2.2) i> l/2+ {x)=e-^°^Y. 



2000 , \ 



e 

71=1 



15 

n 2 



Here there is more evidence that the amplitude of the wave at p = \ is decreasing 
with x — log(fc). The experiments of Fig. 3 give in any cases a stronger evidence: 
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Figure 3. Plot of ip p for p = 1, |, §, |, |, |, ^ up to x = 30, in 
order of increasing amplitudes 




Figure 4. Plot of 4>i/ 2 + 



that for p > | the amplitudes of the waves are decaying, and thus are bounded in 
amplitude by a constant. This is a symptom of the absence of nontrivial zeros in 
the critical segment | < p < 1. 

In the last experiment we set p = | and compare ^3/4 with the asymptotic 
expression of Baez-Duarte: for the trivial zeros we set p = | in the above formula, 
for the nontrivial zeros (the two we consider) we keep the same value of r 3(z) but we 
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assume that their real part is 9t(z) = |. The plot in Fig. 5 of the function ^3/4 fa) 
and of 53/4(2;) + n fa) + r 2 fa) are clearly different: in V'3/4 there is the trace via the 
Mobius function of where the nontrivial zeros are located and thus the amplitude 
is decaying. In the second function, the two considered zeros are supposed to have 
9tfa) = § and the wave which appears seems to have a constant amplitude as in 
the case V'1/2 which of course would be sufficient to ensure the truth of the RH. 




Figure 5. Plots of the functions ^3/4 (x) [red] and 33/4 (x)+rifa) + 
7*2 (x) [green] 

In Appendix we analyse a (weak) stability property of our results obtained with 
N = 2000 in the Mobius function and give some indications why the waves for 
p = I should be decaying, thus ensuring more credibility on the absence of zeros of 
the Riemann Zeta function in the segment | < p < 1. 

3. Conclusions 

In this work we have analyzed numerically the behaviour of the Riesz and Hardy- 
Littlewood wave (the critical function) in some details for the case a = ^ an d = 4 
in the region up to about k = lO'OOO milliard for various values of p in the critical 
segment h < p < 1. In the variable x = log(fc) up to 14 oscillations have been 
detected whose amplitude has been compared with the one calculated with the 
expansion of Baez-Duarte using the trivial zeros and only the two lower nontrivial 
zeros. The agreement is satisfactory and the results give some indication in the 
direction to believe that at least for SH(s) = p > | there are no nontrivial zeros 
of C(s) since in the representation seems to remain bounded. In addition, we 

have given some evidence that a slow decay of Ck like l ° 6 [ k ^ in between to the decay 

-V (RH for the model) and -\ (RH for the model with simple zeros) is possible in 

k 4 k 4 
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the range log(fc) < 30. A further study in this direction but by means of two new 
representations of the Zeta function with coefficients bk and dk and other type of 
oscillations will be presented in the near future. 



Appendix A. 

We consider the critical function ^3/4 ( x ) obtained with TV = 2000 (maximum 
argument in the Mobius function appearing in the Baez-Duarte definition of the 
Cfe). We will suppose that the numerical results are given with good accuracy; we 
now ask: if we increase TV from 2000 up to 10 6 in a ideal experiment, what will be 
the change of the critical function in the range x < 30? 

2000 , , 

^3/ 4 (x; N = 2000) = e% x V ^Je"S 

n=l 712 
10 6 

Vv 4 (*;A^10 6 ) = e^£4re-S. 

n=i n2 

The difference A between the two functions is bounded (1^(^)1 < 1) by: 

io 6 1 
A < e^ x } 

i=200i 

If we ask that A will be smaller then say 10~ 6 time 0.015 which is about the 
value of the amplitude of the wave in the range x < 30, obtained with N = 2000, 
we have: 

x -jr 15 
A < e^ x e~^{C,{ — ) -({ — ;N = 2000)) < 0.015 • 10~ 6 . 

The difference between the Zetas is estimated to: 

oo 

1 , 2 13 2 



/ -^-d.: ----- ---72( 100 - —10 

J X 2 

2000 



-26 



m ' 13 ~~~ 65' 



And the inequality takes the form: 
27 x 2 

—x - el^ + log(— ) - 26 log(10) + 6 log(10) - log(0.015) < 0, 
16 65 

with the solution x < 27. Thus for x < 27, the amplitudes will change at most 10~ 6 
time of its value 0.015. This shows some stability in the numerical experiments as 
TV increases in a ideal experiment. Of course this is independent of how many zeros 
are employed in the Baez-Duarte estimation. 

Finally, application of the crude inequality as above shows [5] that for larger 
values of x, with N sufficiently big, the numerical results of an ideal experiment 
using the Mobius function in Ck will give a value of the amplitude for example 
smaller then \ of 0.015 (still with N = 2000 at x > 35, with N = 10 9 at x > 88), 
indicating that the value of the amplitude becomes possibly smaller. This is an 
indication in the direction to believe that for p > | there are no nontrivial zeros of 
the Zeta function. 
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To end up with the Appendix it should be remarked that in our experiments we 
have used the formula for the ct [4], given by: 



instead of the formula: 



u(n) 

c fe = > ^U. e ~^. 

n=l 



71=1 V 

Again, as above, the crude inequality \n{n)\ < 1 may be used to show that the 
difference between the two functions i.e. the fluctuations become smaller as k get 
bigger and depends on a and (3. In fact they behave unconditionally as: 



If 1 \ 



i \ fe c 



< 



n=l 

To see that, let A = \dk — Ck\ then: 



^ — Z_v r, a I 6 ( ^ r, 3 ) ) — Z_v # 



n— 1 \ x ' / n— 1 

since e ^ > (1 ^) and passing to the continuous variable x, the contribution 

of the second integral is given by [4]: 

if, 1 \ \ ir(y)r(Hi) 



— 1 77 dx = - 



/3r(^i + fc + i)' 



while the first is given by: 



oo k 

~ e ^ 1 a — 1 

- — dx = —L—n^-r- -). 

f3k^ P 



OO 

/ 



1 

At large k the fluctuation behaves like the difference, i.e. as: 

C 



A < 



For the model under consideration the decay is as -§r and is stronger then in 

fe"S" 

the usual Riesz case (a = j3 = 2) where an early more detailed calculation gives a 
decay like [6]. 

fe2 

Finally it should be added that the general upper bound for A is related to the 
discrete derivative of the Baez-Duarte coefficients given by: 

'" M(«) (<\ 1 \ k (\ 1 X ' 



c k (a + f3,f3), 



n=l \ v 7 v 

n=l x ' x 

which unconditionally are bounded by as above [4]. 

k 3 
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In the same way 



d_ 

dk 



n— 1 n=l 




n a+l3 



k 



which gives the same decay since the function is equal to Ck{a + (3, (3) as above. 
At large k we also have [4] : 



[1] S. Albeverio, C. Cebulla, Miintz formula and tero free regions for the Riemann Zeta function, 
2006, Bull. Sci. Math. 

[2] L. Baez-Duarte, A new necessary and sufficient condition for the Riemann Hypothesis, 2003, 
math.NT/0307215 

[3] L. Baez-Duarte, A sequential Riesz-like criterion for the Riemann Hypothesis, 2005, Inter- 
national Journal of Mathematical Sciences, pp. 3527-3537 

[4] S. Beltraminelli, D. Merlini, The criteria of Riesz, Hardy-Littlewood et al. for the Riemann 
Hypothesis revisited using similar functions, 2006, math. NT/0601138 

[5] S. Beltraminelli, D. Merlini, Riemann Hypothesis: The Riesz-Hardy-Littlewood wave in the 
long wavelength region, 2006, math. NT/0605565 

[6] J. Cislo, M. Wolf, Equivalence of Riesz and Baez-Duarte criterion for the Riemann Hypoth- 
esis, 2006, math.NT/0607782 

[7] M. Coffey, On the coefficients of the Baez-Duarte criterion for the Riemann Hypothesis and 
their extensions, 2006, math-ph/0608050 

[8] K. Maslanka, Baez-Duarte criterion for the Riemann Hypothesis and Rice's Integrals, 2006, 
math.NT/0603713 

[9] M. Wolf, Evidence in favor of the Baez-Duarte criterion for the Riemann Hypothesis, 2006, 
math.NT/0605485 

S. Beltraminelli, CERFIM, Research Center for Mathematics and Physics, PO 
Box 1132, 6600 Locarno, Switzerland 

E-mail address: stef ano .beltraminelli@ti . ch 




a Poisson like distribution for the coefficients Cfe. 
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